Particle self-bunching in the Schwinger effect in spacetime-dependent electric fields 
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Non-perturbative electron-positron pair creation (Schwinger effect) is studied based on the DHW 
formalism in 1+1 dimensions. An ab initio calculation of the Schwinger effect in the presence of 
a simple space- and time-dependent electric field pulse is performed for the first time, allowing for 
the calculation of the time evolution of observable quantities such as the charge density, the particle 
number density or the total number of created particles. We predict a new self-bunching effect of 
charges in phase space due to the spatial and temporal structure of the pulse. 
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Introduction: The vacuum of quantum electrody- 
namics (QED) is unstable against the formation of many- 
body states in the presence of an external electric field, 
manifesting itself as the creation of electron-positron 
pairs [IHll • This effect has been a long-standing but still 
unobserved prediction as the generation of near-critical 
field strengths ^ 10^*V/m has not been feasible so 
far. Due to the advent of a new generation of high- 
intensity laser systems such as the European XFEL or 
the Extreme Light Infrastructure (ELI), this effect might 
eventually become observable within the next decades. 

Previous investigations of the Schwinger effect in the 
presence of time-dependent electric fields ji-14|, space- 
dependent electric fields 3; 3 as well as collinear elec- 
tric and magnetic fields |20l - l22| led to a good understand- 
ing of the general mechanisms behind the pair creation 
process by now. However, realistic fields of upcoming 
high-intensity laser experiments showing both spatial and 
temporal variations have not been fully considered yet. 
Only recently it became possible to study the Schwinger 
effect in such realistic electric fields owing to recent theo- 
retical progress as well as due to the rapid development of 
computer technology. Specifically, the Dirac-Heisenberg- 
Wigner phase-space formulation of QED in the presence 
of an external electric field (23l - [26l (DHW formalism) has 



attracted interest again It provides a real-time 

non-equilibrium formulation of the quantum production 
process. Also, a one-to-one mapping between the DHW 
function (phase-space formalism) and the one-particle 
distribution function (quantum kinetic formalism) exists 
in the limit of a spatially homogeneous, time-dependent 
electric field. 

The Schwinger effect in the presence of an arbitrary 
spacetime-dependent electric field is properly described 
by the DHW formalism in the form of a partial differ- 
ential equation (PDE) system for the irreducible compo- 
nents of the DHW function. The numerical solution of 
the PDE system allows for the calculation of any observ- 
able quantity in terms of the irreducible components. 

In the present work, we consider a simple model for 
a sub-attosecond high-intensity laser pulse in standing 



wave mode with finite extension. In the focus of the 
beam, pair production along the direction of the electric 
field gives the dominant contribution to the Schwinger ef- 
fect. Ignoring particle momenta orthogonal to this dom- 
inant direction, the system reduces to a 1+1 dimensional 
setting, which is studied for the first time here and solved 
numerically (30| . 

Formalism: Following the fundamental work of (23 |. 
we start with the gauge- invariant equal-time commutator 
of two Dirac field operators: 

<Pix,y,t) ■.= U{x,y)[^ix-y/2,t),^ix + y/2,t)] , (1) 

with X denoting the center-of-mass and y the relative 
coordinate. Here, the Wilson-line factor which ensures 
gauge invariance is chosen along a straight line: 

Uix, y) - cxp (^le J ^^^^ d^A^x + ^y, t) . (2) 

The vector potential A{x,t) is treated as classical mean 
field, i.e. photon fluctuations are neglected. This approx- 
imation is well justified for the pair-production process 
in QED. Tree-level radiation reactions which might play 
a sizable role for strong fields according to recent inves- 
tigations (siMsst are also neglected in this work. 

Taking the vacuum expectation value (Q\ ^(x, y, t) |f2), 
we trade y for a kinetic momentum variable p by a Fourier 
transformation. This defines the DHW function: 



dye-^py{nmx,y,t)\n). (3) 



Due to the fact that >V(a;,p, t) is in the Dirac algebra, it 
may be decomposed iir terms of its Dirac biliirears: 



W{x,p,t) = i[s + z7^p + 7''^ 



(4) 



with irreducible conrponents transforming as scalar 
s{x,p,t), pseudoscalar p{x,p,t) and vector v^{x,p,t). 
For brevity, these components will later on collectively be 
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denoted as w{x,p,t). The derivation of the correspond- 
ing equations of motion follows that in 3+1 dimensions 
24 130| and yields the following hyperbolic PDE system: 



2pp = 



1:^ 



■ A] p + 2p s 

with the pseudo-differential operator 

nl/2 



, 
, 
-2mp , 
2mv , 



e I d£_E{x 

1-1/2 



^6 



(5) 
(6) 
(7) 
(8) 



(9) 



Along with uj{p) — \/ vn? + p^, the appropriate vacuum 
initial conditions at asymptotic times t^ac — ^ —oo are 



uj(p) 



and Vvac(p) 



(10) 



The irreducible components are not directly observable, 
however, they constitute the observable quantities which 
can be derived from Noether's theorem. For our purpose, 
the charge Q{t) as well as the energy of the Dirac particles 
£{t) are of special interest: 



Q(i) =ey drvo(a;,p,<) , (H) 
£it) = dr[ms{x,p,t) +pw{x,p,t)] , (12) 



with dr = da;dp/(27r) denoting the phase space vol- 
ume element. The integrands q{x,p,t) = vo{x,p,t) and 
e{x,p,t) = [ms{x,p,t) + pv{x,p,t)] can be viewed as 
pseudo-charge density and pseudo-energy density, respec- 
tively. Due to the fact that we are considering a quantum 
theory, it is more appropriate to consider the momentum 
space marginal distributions: 



e{p,t) 



dx 
dx 



q{x,p,t) 



(13) 



[■ms{x,p,t) + pv{x,p,t)] . (14) 



Requiring that the total energy of the Dirac particles 
should be calculable by integrating a particle number 
pseudo-distribution n(x,p^t) times the one-particle en- 
ergy Lo{p), it is also useful to introduce the momentum 
space particle number densities 



n{p,t) 



dx 

W) 



n{x,p,t) , 



(15) 



with 



n{x,p, t) 



m[s{x,p,t)- 

85 vac (p)] +P['^{x,p,t)- ip)] 



(16) 



The vacuum subtractions account for a normalization of 
the density relative to the vacuum Dirac sea. Accord- 
ingly, the total number of created particles reads: 



(17) 



N{t) = / dp nip, t) 



The PDE system Eqs. © - calls for further rewrit- 
ings or even approximations as arbitrarily high momen- 
tum derivatives have to be taken into account in general: 

(a) Full solution in conjugate space: As the momentum 
p appears linearly in the PDE system Eqs. ([5|) - (jS]), we 
can transform these equations to conjugate y space. As 
a consequence, A{x,p,t) transforms into a function of y 
as well: 

f ^eWA(x,p,t) = -ley f'^ d^Eix + ^y^t) , 

(18) 

resulting in an exact, first order PDE system. 

(b) Leading order derivative expansion: The simplest 
approximation is to expand A(a;,p, t) in a series with re- 
spect to the spatial variable. Requiring that [l^]: 



\E{x,t) 



dm{x,p,t)] 



> 



d^wix,p,t)] I 



j.\E"ix,t)^ 



(19) 



it is well justified to neglect the higher derivatives: 

Aix,p,t)^eEix,t)-§^ , (20) 

yielding an approximate, first order PDE system. 

(c) Local density approximation: Approximations can 
also be constructed on the level of the marginal distribu- 
tion n{p,t). Given an electric field E{x,t) = Eog{x)h{t), 
and assuming that the spatial variation scale is much 
larger than the Compton wavelength A ^ Ac, it is well 
justified to locally describe the Schwinger effect at any 
point X independently. We then define the particle num- 
ber quasi-distribution in local density approximation as: 



nioc{x,p,t) := 2T{p,t;x) 



(21) 



Tij), t; x) denotes the one-particle distribution function 
which is found by solving the quantum Vlasov equation 
(33 . [35I at any fixed point xaxcd for a time-dependent 
electric field E{t) = Eog{xtixcd)h{t). Accordingly: 



niocip,t) 



dx 



n\oc{x,p,t) 



(22) 



Results: Our idealized model for a spatially and tem- 
porally well-localized laser pulse in a standing wave mode 
is parameterized by the electric field: 



E{x, t) = Eq exp I 



2A2 



sech' 



2 / 1 ^ 



(23) 



with T and A denoting the characteristic time and length 
scale, respectively. We choose the parameters t — 10/to, 
Eq = O.S-Ecr in this investigation, corresponding to an 
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intense sub-attosecond pulse. As the spatial extent as 
well as the total energy of the electric field of the pulse 
decrease with A, if all other parameters are held fixed, it 
is convenient to disentangle this trivial scaling effect and 
investigate scaled quantities for better comparability: 

and A^W :^ ^ . (24) 

Full solution vs. approximations: In Fig. [1] we compare 
the asymptotic value n{p, t — oo) of the full solution with 
the leading order derivative expansion as well as with the 
local density approximation for different values of A. The 
difference between the various results is rather small for 
broader pulses. As the various approximations are in 
good agreement with the full solution, the pair creation 
process can indeed be considered as taking place at any 
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FIG. 1. Comparison of n{p, t oo) for the full solution (solid) 
with the l.o. derivative expansion (dashed) and the local den- 
sity approximation (dotted) for r = 10/m, Eq = 0.5i5cr- 



point X independently in this regime. For decreasing A, 
however, the various results differ substantially. 

As expected, the leading order derivative expansion 
becomes worse for small A. Whereas a previous study 
of higher derivative terms signalled a potential failure 
at large momenta we here observe a breakdown of 
this approximation for small momenta p/m — > 0. For 
larger A, the dominant momenta are still well approxi- 
mated, but for A approaching Ac, the truncation arte- 
facts overwhelm the physical values. Also the fact that 
the particle density n(p, t — >■ oo) acquires negative values 
in the derivate expansion signals a clear breakdown of 
this approximation for small momenta. The local den- 
sity approximation fails in a different respect: The peak 
momentum of the full solution is shifted to smaller val- 
ues for decreasing A which is not reflected by the local 
density approximation. 

Particle number density: In Fig. [21 we investigate the 
behavior of the full solution n{p,t — > oo) for different 
values of A. A decreasing A involves a shift of the peak 
momentum to a smaller value: The value of the acceler- 
ation by the electric field depends on the actual position 
such that the field excitations feel a varying acceleration 
when moving through the electric field. Accordingly, the 
field excitations are less accelerated for narrow pulses. 

Morover, the shape of n(p, t — oo) becomes higher and 
narrower for decreasing A, at least for A > 4 Ac. This is a 
self-bunching effect caused by the spatial inhomogeneity: 
Excitations which are created with high momenta are 
accelerated for a shorter period as they leave the field 
rapidly. By contrast, excitations which are created with 
small momenta stay longer inside the field and are ac- 
celerated for a longer period. Accordingly, the created 
particles are bunched into a smaller phase space volume. 

For A < 4Ac, however, the height of n{p,t — s- oo) de- 
creases again as more and more field excitations gain too 
little energy in order to finally turn into real particles. 
For A = Ac, the energy content of the electric field is 
ultimately so small that none of the vacuum fluctuations 
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FIG. 2. Comparison of n(p,t oo) for r = 10/m, Eo = 
0.5iJcr and different values of A. Note that n{p, t — >■ oo) = 
for A = Ac- 
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FIG. 3. Comparison of N{t — )■ oo) for the full solution (solid) 
with the l.o. derivative expansion (dashed) and the local den- 
sity approximation (dotted) for r = 10/m, Eq = 0.5-Ecr- 



These results suggest further studies of the Schwinger 
effect in realistic space- and time-dependent electric fields 
in 3+1 dimensions. The goal is to consistently describe 
the Schwinger effect beyond the mean field level by tak- 
ing into account photon corrections to the background 
electric field and subsequent collision and radiation pro- 
cesses. In the long run, we expect the self-bunching effect 
to play an important role in the generation of taylored 
electron/positron beams. 
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eventually turns into real particles. This observation is in 
good agreement with previous studies of space-dependent 
electric fields E{x) |l5l4l7j : The pair creation process is 
expected to terminate once the work done by the elec- 
tric field over its spatial extent is smaller than twice the 
electron mass. As the pair creation process occurs at 
time scales of the order of the Compton time tc = l/m, 
which is smaller than the time scale of the electric pulse 
r = 10/m, this estimate should be reasonable in our case 
as well. The corresponding estimate for the pair creation 
process to terminate for Eq ~ 0.5Ecr is in fact in good 
agreement with our results: 



A<^J-Ac^l.6Ac. 
Eo \ n 



(25) 



Number of created particles: In Fig. [3] we compare the 
asymptotic value N{t — > oo) obtained from the full solu- 
tion with the leading order derivative expansion as well 
as with the local density approximation for different val- 
ues of A. Again, we observe good agreement between the 
full solution and the various approximations for large A, 
however, substantial deviation for small A. Most notably, 
only the full solution shows the sharp drop of N{t oo) 
for small A in accordance with Eq. ([55]) . 

Conclusions: We have presented an ab initio real- 
time calculation of the Schwinger effect in the presence 
of a simple space- and time-dependent electric field pulse 
in 1+1 dimensions, showing various remarkable features: 

Most notably, we observe a new self-bunching effect in 
phase space which can naturally be interpreted in terms 
of the space and time evolution of the quantum excita- 
tions. The pair creation process eventually terminates 
for spatially small pulses once the work done by the elec- 
tric field is too small in order to provide the rest mass of 
an electron-positron pair. Whereas the derivative expan- 
sion is quantitatively able to signal these self-bunching 
effects, the local density approximation fails to describe 
these important properties. 
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